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Radial equation for spherically symmetric potential 



The SE in 3D in spherical coordinates is 



(22-1) 



using the ansatz if)(r) = R(r)Y(9, 0), and inserting for the angular function an eigen- 
function 

= Y lm {9,<j>) = (8,<f>\l,m), (22-2) 

we have, using \?Yi m {9 , 0) = % 2 l(l + \)Yi m {9, <f) after dividing by Y\ m for the radial 
equation, 



2m \dr 2 ' r dr 



2 2 9 
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k 2 l(l + 1) 
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+ V(r) 



i? n ,(r) = E nl R n i(r) 



(22-3) 



Here, we have added two subscripts n, I to the radial wavefunction R(r) and the 
eigenenergy E because the SE for the radial part of the wavefunction depends on the 
total angular momentum / of the 3D wavefunction ip(r). 

Note. The z-component of angular momentum L z , and the corresponding magnetic 
quantum number to, do not appear in the radial equation. 

We can define an /-dependent effective potential, 

h 2 i{i + 1) 



Kff,i = V{r) + 



(22-4) 



2mr 2 

where the additional term is the centrifugal barrier for a particle with angular mo- 
mentum 

(L 2 ) = Ti 2 l(l + 1). (22-5) 

The radial equation can be brought into a more familiar-looking form by introducing 
a new function: 

u(r) = rR(r) 
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u{r) 



Then, 
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and the radial equation is 



n 2 1 d 2 u 



2m r dr 2 



+ 



h 2 l(l + 1) 
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+ V(r) 



u(r) Tp u ( r ) 
= hi 



or 



h 2 d 2 h 2 i{i + \) _ _ / / 



2m dr 2 
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u n i(r) 



= E nl u nt {r) 



(22-11) 



(22-12) 



This equation for u(r) = rR(r) has the same form as the ID SE in the effective 
potential 

h 2 i(i + i) 



Wr) = V(r) + 



(22-13) 



2mr 2 ' 

but with slightly different boundary conditions. Therefore, u(r) looks like an anti- 




Figure I: u(r) = rR{r) has the same form as the ID SE in the effective potential 
Kff,«( r )> but with slightly different boundary conditions. 

symmetric solution in all space. Consequences are, e.g., that since an antisymmetric 



M 




Figure II: u(r) looks like an antisymmetric solution in all space. 



bound state does not always exist in ID, that a bound state does not always exist 
in 3D (in contrast to ID, where a symmetric bound state always exist in a potential 
well). 3D wavefunctions u(r) are like antisymmetric ID wavefunctions in the effective 
potential 

Vl {r) = V(r) + ^±11. (22-14) 
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Hydrogen atom 



Ze 2 

V(r) = — - — > (and the radial equation is) (22-15) 

UQ2 d 2 Ze 2 k 2 l(l 



+ — - — - E )u(r) = (22-16) 
2m or 1 Aner 2mr l J 

We introduce a dimensionless position coordinate p by p 2 = 8n jj,f\ r 2 , and define for 
E < 




Ze 2 8m h 2 



16nehr \ \E\ \ 8m\E\ 



Ze I m 1 



A 



(22-17) 

(22-18) 
(22-19) 



The equation can be written as 



&_ fX 1 1(1 + 1) 

dp 2 \p 4 p : 



2 „+ ----^_^)„ = (22-20) 



with p = sj^r, A = S^/^ = ^a^/fj, where a = ^ » ^ is the dimen- 

sionless fine structure constant. To solve this equation, we proceed as for the HO: We 
write a Taylor-expansion solution after having factored out the correct asymptotic 
behavior. 



For very large p we have 



For very small p, 



d 2 1 

— u = —u (22-21) 
dp 2 A 

u(p) oc e~* p (22-22) 



£«=^« (22-23) 
dp 2 p 2 

u(p) oc p l+1 (22-24) 



Consequently, we try a solution of the form 



u(p) = s(p)p l+1 e-^ p (22-25) 
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u 
u 



'(p) = (s'(p)p l+1 + s(p)(l + l)p l - \sp l+1 )e-^ 
"(p) = [ S 'V +1 + 2(1 + l)s'p l + s(l + 1) Ip 1 - 1 
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(22-26) 
(22-27) 
(22-28) 

(22-29) 

l + l 1 " 

s + -s 

p 4 

(22-30) 
(22-31) 



Inserting this into (??) leads to 



s" + 



2(1 + 1) 



l)s' + 



([^_[±1 + 1 A l_fcl) )s = (22 . 32) 
p 2 p 4 p 4 p 2 



2/ + 2 
P 



. A — Z — 1 

s' + s = 

P 



To solve this differential equation, we write a Taylor expansion about p = 0: 



oo 

>(p) = 5^ fl fep A 

A:=0 

oo 



Jfc-2 



^a k k(k-l)p k 

k=0 

oo 

Y j a k+ 2(k + 2)(k + l)p k 
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(22-33) 

(22-34) 
(22-35) 
(22-36) 
(22-37) 
(22-38) 

(22-39) 



(21 + 2)J2 a k+ 2(k + 2)p k - ak + i(k + l)p k 

k=0 

oo 

s = (A - I - 1) V a k+1 p k 

P — 

F fc=0 

which substituted into (??) results in 
Y,Pk{(k + 2)(k + l)a k+2 + 2(l + l)(k + 2)a k+2 + (\ — l — l — k— l)a k+1 j = (22-40) 
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This must vanish term by term, so we obtain a recursion relation 
(k + 2){k + 21 + 3)a k+2 = (k + l + 2 - \)a k+1 



(22-41) 



or 





recursion 


ctfc+i k + I + 1 - A 


relation for 


a k (k + 1) (k + 2(1 + 1)) 


expansion 


coefficients 



(22-42) 



If the series does not break off somewhere, we will have for large k, a k oc \o-k~i or 
a k oc |p which gives a growth s(p) oc e +p , which is not acceptable for u(p) = s(p)e~2. 
Consequently, we require the series to terminate, which implies A = k + 1 + 1 for some 
L. Let us call n r = k the integer with that property. It is cutomary to define the 
principal quantum number as 



n = n r + I + 1 , 
where n r > 0, so n > 0, so n > I + 1, n integer, and 

Ze 2 I m~ 

An 



(22-43) 



AuehM 2\E„ 



2\E„ 



n 



(22-44) 

(22-45) 
(22-46) 



Consequently, the eigenenergies of the hydrogen atom are 



eigenenergies of 
hydrogenlike atoms 



(22-47) 



This is the same energy eigenspectrum as obtained from the Bohr formula. 
Note. There are important differences: 

• The principal quantum number n = n r + I + 1 is really the sum of the radial 
quantum number n r and the total angular momentum quantum number /. 

• We have obtained the full radial and angular distribution of the electron, which 
generalizes the classical concept of an orbit. 
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First few radial functions 



with the Bohr radius 
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(22-48) 
(22-49) 
(22-50) 
(22-51) 
(22-52) 



(22-53) 



1. n r = / = 0, n = m = A, a\ = 



1 / Zr\ Za 

u(r) = Cpe 2P = Ci( — )e a o 
R(r) = ^= C 2 e~^ 



(22-54) 
(22-55) 



Note. The probability to find the electron between r and r + dr is given by 
r 2 \R(r)\ 2 dr = \u(r)\ 2 dr. 



2. (a) n r = 1, I = 0, n = 2 = A 



Oi 

a 
«2o(^) 



1-2 



1 
2 



V 2^ a V 2a ; 



(22-56) 
(22-57) 
(22-58) 
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(b) n r = 0, I = 1, n = 2 = A 



Oi 

a 



-> a : = (22-59) 

(22-60) 



U21 (r) = CpV^ = C"^^ e"^ 



i22i(r) = C"f— (22-61) 

i?2o — R n =2,i=o and i?2i = R n =2,i=i are different states that have the same 
eigenenergy. The occurrence of different eigenstates with the same energy, 
(or in general quantum number) is called degeneracy. 
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